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ABSTRACT 

The recent discoveries of ultra-faint dwarf (UFD) galaxies in the vicinity of the 
Magellanic system supports the expectation from cosmological models that such faint 
objects exist and are numerous. By developing a mass model of the Local Group and 
backwards integrating the Magellanic Clouds’ present kinematics, we find that the lo¬ 
cations of these UFDs are consistent with those predicted if previously associated with 
the Large MC as part of a loose association. We further demonstrate how these satel¬ 
lites are likely to have been processed by the Galactic hot halo upon accretion, with 
the implication that ongoing detections of extremely gas-rich objects on the periphery 
of the Galaxy and without clear stellar counterparts are analogous to the progenitors 
of the gas-deficient UFDs. Our model allows us predict the locations of other puta¬ 
tive Magellanic satellites, and propose how their distribution/kinematics provide a 
novel constraint on the dynamical properties of the Galaxy. We also predict that the 
stripped metal-poor H I, previously associated with these UFDs, lies coincident with 
but distinguishable from the extensive Magellanic Stream. 
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1 INTRODUCTION 

The last decade has seen the proliferation of resolved or 
candidate satellites in the Local Group (LG), whether 
by wide-field imaging surveys and/or blind HI surveys 
(e.g Simon & Geha 2007; McConnachie 2012; Adams et al. 
2015), in apparent fulfillment of the prediction for adun- 
dant dark substructure by the standard ACDM cosmo¬ 
logical model (Klypin et al. 1999). Many constitute ultra 
faint dwarf (UFD) galaxies with, in some cases, as little as 
10 3_4 Mq in stellar mass and mass to light ratios of >1000 
Mq/Lq. Ongoing and future observations of such systems 
provide great insight as to the baryonic processes dictating 
the faint end of the luminosity function (Brown et al. 2015; 
Meyer et al. 2015). 

These same observations can resolve another tension 
with the ACDM model, insofar as the spatial distribution 
of most LG satellites is not isotropic as expected but ex¬ 
hibits a correlation in phase-space (Kroupa, Theis & Boily 
2005); a notable recent interpretation of which is that these 
satellites are the tidal remnants of a past M31 merger event 
(Hammer et al. 2013). 
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An alternative explanation for this coherence lies 
in a group infall scenario, centred on the Large Mag¬ 
ellanic Cloud (LMC; Lynden-Bell 1976; Bekki 2008; 
D’Onghia & Lake 2008). Cosmological simulations of the 
local volume suggest Galactic-type halos preferentially ac¬ 
creted satellites from 10 10 ~ 12 Mq subgroups 5 to 8 Gyr ago 
(Wetzel, Deason & Garrison-Kimmel 2014), in which case 
the uncharacteristically blue and gas-rich Large MC is a 
prime candidate for a presumed subgroup central galaxy 
(Tollerud et al. 2011). Studies of known dwarf associations 
on the periphery of the LG (Tully et al. 2006) suggest how¬ 
ever an incompatibility between their distribution and that 
of the majority of LG satellites (Metz et al. 2009). 

Recent independent detections of stellar overdensi¬ 
ties in the DES survey by Koposov et al. (2015, K15), 
Bechtol et al. (2015, B15) and Kim & Jerjen (2015) have 
yielded the best candiates for Magellanic satellites thus far, 
with magnitudes —6.6 < My < —2.6 and lying within a 
relatively small footprint close to the LMC. It is as yet un¬ 
certain if they fulfil an expectation from cosmological sim¬ 
ulations that LMC-analog Galactic satellites would be ac¬ 
companied by fainter objects lying in its vicinity of the 
sky (Sales et al. 2011). By identifying LMC-mass galaxies 
accreted onto MW-type halos in the ELVIS suite of N- 
body simulations, Deason et al. (2015) suggest that surviv- 
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Figure 1 . Galacto-centric radius of the LMC for orbits BOUND 
and INFALL, whose initial conditions differ only by the assumed 
present-day MW mass/c; the respective 1 <t spread in UFD model 
orbits with respect to the LMC are shown. The present radii of 
the observed UFDs (K15/B15) are shown with black dashes where 
the top-most corresponds to Eridanus 2; (Inset) LMC orbits in 
x-y space (coplanar with and centred on the MW disc). 

ing LMC-type satellites account for as much as 7 percent of 
the LG inventory, although the specific association of these 
UFDs and others (e.g Leo II/IV/V) hinges on whether the 
LMC has just passed its first pericentre. 

In this paper, we establish the likelihood that these 
UFDs were previously associated with the LMC as part of 
a loose association using the latest estimates for the dy¬ 
namical/kinematic properties of the Galaxy/LMC, which af¬ 
fords the opportunity to consider group processing specific 
to the local environment. The absence of H I-counterparts 
to these UFDs (Westmeier et al. 2015) leads to demonstrate, 
therefore, an evolutionary link between these stellar systems 
and the so-called gas-rich dark galaxy candidates, such as 
recently-discovered Leo P and AGC198606 located on the 
periphery of the Galactic halo (Adams et al. 2015). 

2 ORBITAL HISTORY OF A HYPOTHETICAL 

UFD PROGENITOR 

2.1 Method 

We adopt a method that builds upon that of Nichols et al. 
(2011), in which a long-term MC-UFD association is de¬ 
duced from tracing the trajectories of the LMC and small 
MC (SMC), with respect to the MW, back by 9 Gyr 
(at which point infall is deemed to have commenced; 
Wetzel, Deason & Garrison-Kimmel 2014), then integrating 
forwards 10 s UFDs models according to a Monte-Carlo tech¬ 
nique, wherein their initial location is randomly sampled 
from an assumed Magellanic Group that is analogous to 
known dwarf associations (Tully et al. 2006). All orbits are 
integrated with the symplectic leap-frog algorithm. 

Using the HST third-epoch proper motions for the MCs 
(Kallivayalil et al. 2013, K13), we trace the past MC tra¬ 
jectories with a backwards-integration method composed 
of potential models of the LG (Murai & Fujimoto 1980), 


as adopted previously to establish the formation of the 
Magellanic Stream (MS; Yozin & Bekki 2014). Our model 
therefore consists of a MW comprising an adiabatically- 
contracted NFW dark halo with dynamical mass Mmw, a 
Miyamoto-Nagai stellar disc with mass Md and a Hernquist 
bulge with mass Mb. The dynamical friction imposed on 
the orbits of the MCs while travelling through the MW 
halo is modelled according to the method established by 
Zenter & Bullock (2003). Our models of the MCs each con¬ 
sist of NFW halos with mass Mlmc (Msmc) and exponential 
discs with mass M d ,LMC (M diSM c)- 

The redshift evolution of the respective NFW 
mass and concentration (c) of the MW is in¬ 
corporated here using merger rate-* relations 
from Fakhouri, Ma & Boylan-Kolchin (2010) and 
Munoz-Cuartas et al. (2011) and assuming the MW halo’s 
virial radii (r v ; r ) scales with mass Mmw as r v i r ocMV 3 . We 
also assume the stellar components of the MW, M d and Mb, 
have grown linearly by 50 percent since z = 2 (Snaith et al. 
2014). 

The stellar components of the MCs are fixed in mass, 
but their halo masses, Mlmc and Msmc, are assumed to 
lose mass as a function of time due to tidal truncation. 
The halo mass of each MC at 9 Gyr ago is thus pre¬ 
scribed the value that corresponds to their current stel¬ 
lar mass according to halo abundance matching methods 
(Behroozi, Wechsler & Conroy 2013); the mass loss rate 
that links these values with those deduced from present day 
observations (e.g. K13) is assumed constant. 

Our model thus comprises various simplifications, many 
of which we find the results to be largely insensitive to. 
The key exception lies in the MW’s total mass and c, in¬ 
sofar as they strongly influence the number of completed 
MW-MC orbits since infall. In this work, we consider a 
BOUND case corresponding to a high-mass MW (Mmw = 
1.6 x 10 12 Mq (Gnedin et al. 2010) and a INFALL case with a 
compact/low-mass MW (Mmw = 0.8 x 10 12 Mq (Kafle et al. 
2014). While each of these estimates are based on stellar 
halo radial velocities, the respective assumptions for den¬ 
sity/velocity anisotropy profiles vary significantly, leading 
to estimates of one and zero completed MW-MC orbits re¬ 
spectively. 

In the absence of dynamical measurements of the ob¬ 
served UFDs, we infer their total masses from neighbouring 
LG objects in the size-luminosity plane (e.g Segue 2/Leo 
V; Simon & Geha 2007; McConnachie 2012), thereby adopt¬ 
ing a stellar mass M dj uFD of 10 3,5 Mq. Abundance matching 
of subhalos in this faint regime is highly uncertain, relying 
on extrapolaton over 5 orders of magnitude; for our UFDs, 
modelled as NFW halos, we adopt a halo-to-stellar mass 
M h /M d | UFD at 9 Gyr ago of 10 5 (Fig. 2), as obtained from 
Garrison-Kimmel et al. (2014) who fit their relations to the 
more luminous LG satellites. 

The z = 0 properties of this LG model are summarised 
in Table 1. The thick lines in Fig. 1 illustrate the respective 
orbits of the LMC in both BOUND and INFALL scenarios as 
deduced from the backwards integration. The shaded regions 
that accompany them convey the la deviation in galacto- 
centric radii (r) of the 10 5 UFDs as a function of time, whose 
respective orbits commence in the vicinity of the LMC at 9 
Gyr ago. The initial UFD locations at this lookback time 
are randomly sampled assuming they trace the isotropic 
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Figure 2. Neutral hydrogen mass as a function of stel¬ 
lar mass for galaxies from the literature (Cannon et al. 2015; 
Onorbe et al. 2015). Also shown is the halo-stellar mass rela¬ 
tion from Garrison-Kimmel et al. (2014), derived from abundance 
matching and extrapolated to the region suspected to be occupied 
by UFD progenitors (grey shaded region). 

NFW distribution of DM in groups, with concentration ~3 
(Lin, Mohr & Stanford 2004) up to distances of ~250 kpc 
from the central LMC (Tully et al. 2006), and with veloci¬ 
ties constrained by an Eddington distribution function. 

2.2 Results 

Fig. 1 illustrates how UFDs in a BOUND scenario are far 
more widely distributed in r than in an INFALL scenario 
on account of significant scattering by Galactic tides upon 
first pericentre ^4 Gyr. This result agrees qualitatively with 
that of Deason et al. (2015); the significant proportion of 
UFD models that remain closely associated with the LMC in 
BOUND is also consistent with prior analytical predictions 
of their dispersion timescale (Bowden, Evans & Belokurov 
2014). 

At present, detections of such faint galaxies will be bi¬ 
ased towards those at low r, at evidenced by the clustering 
of the K15/B15 sample at r<100 kpc. While it is therefore 
difficult to support a particular orbit scenario based on the 
observed UFD locations, Fig. 3 demonstrates our prediction 
for the location for further discoveries, where the modelled 
UFD locations at zero lookback time are compared with 
those lying in the DES footprint, in an equatorial projec¬ 
tion. The probability distribution of UFDs as a function of S 
reveals those in BOUND peak at low and high r, while those 
in INFALL are almost equally distributed along the (more 
than ~100 kpc) span of the Magellanic Stream (Fox et al. 
2014), many of which will likely be obscured by it. 


3 GAS LOSS TO THE HOT GALACTIC 
CORONA 

Faint galaxies observed on the periphery of the MW halo 
and in extragalactic surveys (e.g. ALFALFA) report gas- 
to-stellar ratios of >100 (Huang et al. 2012; Cannon et al. 
2015), and the ongoing UV emission from others is proposed 
as a viable means for future discoveries (Meyer et al. 2015). 
By contrast, the UFDs of K15/B15 (with the possible excep¬ 
tion of the anomalously distant and blue Eridanus 2) appear 
to have no associated H I (Westmeier et al. 2015), and the 



Figure 3. (Left panel) Present locations of the LMC/SMC 
(red/blue circles respectively) and their past orbits (BOUND sce¬ 
nario), projected in equatorial coordinates; also shown is the DES 
footprint (red shaded region), approximate location of the Magel¬ 
lanic Stream (MS; blue shaded region), and the present locations 
of observed UFDs (filled circles) coloured according to their he¬ 
liocentric distance. Similarly, the average location of UFD mod¬ 
els binned in distance increments of 40 kpc (up to 400 kpc) are 
shown with filled squares coloured according to their mean dis¬ 
tance; (right) probability distribution of UFDs as a function of 
declination <5 


Table 1. Summary of a) properties of observed UFDs; b) Initial 
UFD model parameters; c) present day MW and MC properties, 
where Mmw an< i c f° r the MW are given for BOUND and INFALL 
(in parentheses). 


(a) UFD Parameters 

Range 

Half-light radius 

0.014 to 0.169 kpc 

My 

-6.6. to -2.7 mag 

Heliocentric distance 

30 to 380 kpc 

(b) UFD Model Initial Parameters 

Value 

Stellar mass M* 

10 3 - 5 M 0 

Gas mass M g 

10 6 M 0 

DM halo mass 

10 8 - 5 M 0 

DM halo radius 

2.5 kpc 

Stellar disc radius 

0.1 kpc 

Gas disc radius 

0.26 kpc 

(c) MW and MC Parameters 

Value 

MW halo mass Mmw 

1.6(0.8)xl0 12 M 0 

MW halo concentration, c 

12 (21) 

MW stellar disc mass 

5.5xl0 10 M 0 

MW stellar disc radius 

3.5 kpc 

MW bulge mass 

lxlO 10 M 0 

MW bulge radius 

0.7 kpc 

LMC (SMC) halo mass 

5xl0 10 (3xl0 9 ) M 0 

LMC (SMC) halo concentration, c 

9 (15) 

LMC (SMC) virial radius 

120 (60) kpc 

LMC (SMC) stellar disc radius 

1.4 (1.1) kpc 
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strong possibility that they have been recently accreted to 
the Galaxy as part of a Magellanic Group would tentatively 
support their recent processing by this environment. 

However, the star formation histories of other Galactic 
satellites suggest that those with the least stellar mass were 
quenched at very high redshift. Brown et al. (2015) combine 
the photometry, spectroscopy and isochrones of six UFD 
Milky Way satellites to find in all cases an ancient stellar 
population, largely formed prior to 2 = 3. Weisz et al. (2015) 
define the epoch of quenching for each of their sample of 38 
satellites as when 90 percent of their respective stars were 
formed, and establish a (possibly underestimated) quenched 
fraction of 0.4 by z = 2, which in many cases predates their 
time of infall to the Galactic system. 

The pre-infall quenching is consistent with the extant 
gas of UFDs escaping their halos upon reionisation of the 
universe, but this scenario contrasts with ultra-high reso¬ 
lution simulations (sufficient to resolve M* ~ 10 4 Mq UFD 
analogues) that have recently probed the effect of the UV 
field on the UFD stellar mass regime and find it to merely 
suppress further gas accretion (Onorbe et al. 2015). The ex¬ 
tended SF histories of several low mass dSphs highlight a 
further feasible scenario (Weisz et al. 2015), wherein early 
quenching by the UV field is followed by gas accretion which 
renews a period of SF that is quenched once again upon in¬ 
fall. 

3.1 Method 

In this section, we investigate the possibility that the ab¬ 
sence of H I associated with the K15/B15 UFDs can be 
explained in terms of an interaction between a previously 
gas-rich disc and the Galactic hot halo alone. We thus pre¬ 
scribe our UFD models with an initial H I-to-stellar mass of 
10 2 ' 5 ; illustrated in Fig. 2, this is motivated by high M*/L* 
systems imaged by the VLA (Cannon et al. 2015) and the 
simulated analogues from Onorbe et al. (2015). 

In an analytical method to establish the gas loss, we 
adopt the same 10 s UFD orbits computed in Section 2, 
which in each case provides their respective galactocentric 
radius r at a given lookback time. The MW’s hot halo den¬ 
sity at r is computed from the modified-/? density model of 
Miller & Bregman (2015); although the normalised hot halo 
mass distribution is fixed, its total mass is assumed to evolve 
in a self-similar manner that follows the evolving Mmw as 
documented in Section 2. 

The gas mass loss rate via ram pressure for each UFD 
orbit is predicated on the condition that: 

pV 2 ^ E(r d )^l| zo , (1) 

where p is the local halo density for a UFD with velocity 
V and a restoring force perpendicular to its disc computed 
from the gas surface density profile E and potential gradient 
Scp/Sz at the in-plane radius rj. Thus we assume the strip¬ 
ping of all H I beyond r<j for which the condition is met. 

A caveat with this method lies in the uncertainties re¬ 
garding the height zo at which the gradient should be com¬ 
puted, and the rate at which stripped gas can fall back upon 
the galaxy. To address this issue quantitatively, we perform a 
single self-consistent N-Body/SPH ram pressure simulation 
for a single UFD orbit (in the BOUND scenario). By calcu¬ 


lating the sum of the galaxy’s SPH gas particles exceeding 
their respective escape speed within the galaxy’s instanta¬ 
neous potential, the corresponding gas stripping rate during 
the critical pericentre passages can be used to calibrate our 
analytical model. 

The details of this numerical apparatus is given in our 
previous studies (Bekki 2014; Yozin & Bekki 2015), but can 
be briefly summarised as follows: a cubic lattice of SPH par¬ 
ticles with periodic boundary conditions are used to model 
the Galactic hot halo local to the galaxy (up to 10 disc 
scalelengths). The self-gravitating galaxy model sits at rest 
within this lattice, while the density and relative velocity of 
the halo SPH particles are varied at each simulation timestep 
(length of order 10 5 yr) according to the aforementioned pre¬ 
defined UFD orbit and the feefa-model density at the corre¬ 
sponding r. 

Although the UFDs are not expected to be rotationally- 
supported, their analytical/simulation models adopt an 
exponential density model for an UFD progenitor with 
a size 0.1 kpc extrapolated from the tight size- 
mass relation (over 10 7 < M*/Mq < 10 11 ) from 

Ichikawa, Kajisawa & Akhlaghi (2012), with a similar pro¬ 
file for the co-planar gas disc, albeit with a factor 2.6 larger 
scalelength (Kravtsov 2013). 

The inclination of these discs with respect to its orbit 
is prescribed to an arbitrary ~37.5 degrees. The analytical 
formula provided above is strictly defined for a face-on (0 
degrees) inclination, although in simulations where the ex¬ 
ternal ram pressure is comparable to the central pressure 
of the galaxy’s gas disc (such as in the scenario considered 
here), Roediger & Bruggen (2006) find the ram pressure effi¬ 
ciency varies only for those inclinations close to edge-on (90 
degrees). By assuming no preference in the infall inclination 
of the alleged Magellanic UFDs, we conservatively account 
for their simulated results by scaling our computed stripping 
rates by a factor 0.7. 

3.2 Results 

The top row of Fig. 4 shows a snapshot from the simulation 
method taken during the first pericentre passage, illustrat¬ 
ing the large H I wake that trails the orbit of the UFD. The 
corresponding gas mass loss is shown in the second row, and 
compared with an equivalent UFD model from our analyt¬ 
ical method and whose properties are calibrated by it; the 
peak H I mass loss rate during the pericentre passages is ~3- 
4x10° M@. The third row compares the gas loss associated 
with ram pressure for a single UFD model constrained to 
the LMC trajectory in both BOUND/INFALL scenarios. To 
account for up to magnitude variations in local halo densi¬ 
ties around simulated stripped satellites (Bahe & McCarthy 
2015), we present for each orbit two runs in which the local 
p is scaled by 1 and 10 respectively. 

This suite of models predict in general a (near) complete 
removal of extant H I by z = 0 for those UFDs proximate 
in location to the LMC, consistent with Westmeier et al. 
(2015). The bottom row of Fig. 4 confirms this result for 
the 10 5 UFD models in a log-scaled histogram of the gas- 
to-stellar ratio-heliocentric distance (r) phase space (where 
p is scaled by random factor /, 0.1 < / < 10). 

The large H I loss rates reported above implies their 
remnants will be similar in both morphology and mass to the 
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^ 10 5 -6 compact High Velocity Clouds that populate the 
MW halo (e.g. Giovanelli et al. 2010). For the instantaneous 
local density and velocity of the UFD, we can calculate the 
destruction timescale of a stripped 10 5 Mq gas blob with 
effective size d = 0.1 kpc and gas density ph ~ 10 19 cm ~ 2 
with the formula of Bland-Hawthorn et al. (2007) i.e. ra — 
ipb/p)°' 5 (d/V). 

If conservatively assuming a velocity V equal to that of 
the UFD, we obtain r of order several Gyr, implying that 
such blobs have feasibly survived to the present day since 
their removal up to 5 Gyr ago, independently of their progen¬ 
itor UFD. While Fig. 3 shows how these blobs most likely 
lie in the vicinity of the MS, it is possible to distinguish 
them from fragments of the MS by their primordial metal 
abundances (e.g Kumari et al. 2015), given that the chem¬ 
ical abdundance and filamentary structure of the Stream 
are consistent with a relatively metal-rich SMC as its origin 
(Fox et al. 2013; Yozin & Bekki 2014). 

The proximity between these the UFD orbits and the 
Stream also implies their mutual liydrodynamical interac¬ 
tion. It is beyond the scope of the methods presented thus far 
to examine this in detail, so we refer to the average gas den¬ 
sity of the Stream (including a substantial ionised compo¬ 
nent), log(nHcm _ 3 )~-1.8 (Fox et al. 2014), and a character¬ 
istic velocity dispersion of group associations (~30 kms^ 1 ; 
Tully et al. 2006) to represent the relative UFD-Stream ve¬ 
locity. The corresponding ram pressure (~pV 2 ) is approxi¬ 
mately two orders of magnitude less than that imposed by 
the Galactic hot halo. Although seemingly insignificant, the 
suspected lifetime of the Stream deduced from tidal mod¬ 
els (e.g. 2 Gyr; Besla et al. 2012) implies a long-term inter¬ 
action that will necessitate an explicit modelling in more 
refined models. 


4 DISCUSSION AND CONCLUSIONS 

In this letter, we have presented a viable formation scenario 
for several of the ultra faint dwarf satellites recently detected 
by Koposov et al. (2015) and (Bechtol et al. 2015). Unlike 
the spatial distribution of dSphs at large in the LG, whose 
origin presently eludes clear explanation within the context 
of ACDM, these UFDs are consistent with their infall as 
part of a loose association centred on the LMC, a scenario 
that agrees well with those of cosmological simulations of the 
local volume (Wetzel, Deason & Garrison-Kimmel 2014). 

If assuming the majority of these UFDs were origi¬ 
nally Magellanic satellites, then our model yields a proba¬ 
bility distribution for the present locations of these UFDs 
(Fig. 3) that predicts in the range of 10 more, most of 
which are lying within the Magellanic Stream. This rep¬ 
resents the second prediction for stellar objects in this re¬ 
gion, where previous attempts to detect the stellar arms 
postulated to have been tidally removed from the MCs 
(Yozin & Bekki 2014) have been thus far inconclusive (e.g 
Ostheimer, Majewski & Kunkel 1997). If such objects were 
found, they could be adopted as valuable distance con¬ 
straints on the Stream itself. 

It is postulated here that overcoming the detection bi¬ 
ases that plague the current K15/B15 UFD sample, and 
revealing the wider population of putative Magellanic satel¬ 
lites, can provide constraints on their combined orbital his- 
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Figure 4. (Top row) Selected snapshots from the N-Body/SPH 
simulation used to calibrate our analytical stripping model, where 
contours represent the gas surface density in log-scaling; (Second 
row) A comparison of the H I-to-stellar mass as a function of 
time for a single UFD model in on a BOUND orbit as computed 
with an analytical method (dashed line), whose parameters are 
calibrated by the results of an equivalent N-Body/SPH simulation 
(solid line); (third row) The same, but for UFDs modelled with 
the analytical method only, and constrained to orbits BOUND 
and INFALL (left and right respectively), where hot halo density 
p at a given r is scaled by 1 and 10; (bottom row) The final 
galacto-centric radius of 10 5 UFDs in both BOUND and INFALL 
scenarios, plotted as a log-scaled histogram (contours represent 
orders of magnitude) in the H I-stellar mass/Heliocentric distance 
phase space. 


tory. By extension, this can provide a novel constraint on the 
wide range of Galactic dynamic masses reported in the lit¬ 
erature; other such constraints have recently been provided 
by models of the tidal formation of the Sagittarius Stream 
(yielding a surprisingly low mass estimate of 4xlO n Mq; 
Gibbons, Belokurov & Evans 2014). 

Our results generally compliment those of Deason et al. 
(2015), who searched for surviving LMC-type satellites and 
their entourage within MW-type hosts of the ELVIS simu¬ 
lation. However, for the range of observationally-motivated 
MW masses considered in this study (0.8 — 1.6 x 10 12 Mq), 
the LMC is an unusual satellite for variety of reasons, not 
least its mass (e.g. Rodriguez-Puebla, Avila-Reese & Drory 
2013). For this reason, the scenario derived from their ELVIS 
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sample is quantitatively different to ours, being biased to¬ 
wards 2 x 10 12 Mq hosts. 

We note however that our non-static MW/MC model, 
like those of Kallivayalil et al. (2013), is subject to a va¬ 
riety of simplifications that require further refinement. As 
noted by Deason et al. (2015), the tidal processing of the 
UFDs prior to and during accretion remains to be consid¬ 
ered, as emphasised not only by the morphological disrup¬ 
tion of other UFDs (Simon & Geha 2007), but also the po¬ 
tentially long-term LMC-SMC binarity as part of a common 
halo (Bekki 2008; Besla et al. 2012). Incidentally, we note 
that although ram pressure acts directly on the UFDs’ gas, 
our single N-Body/SPH simulation demonstrated how the 
significant shift in baryons for such large Mhi/M* (~100) 
drives a morphological change manifesting in the elongation 
of the disc in the direction of its H I-wake. 

We have assumed the effects of reionisation in the UFD 
mass regime are separable to those as part of group process¬ 
ing (Wetzel, Deason & Garrison-Kimmel 2014), and thus do 
not preclude large gas masses upon infall (Onorbe et al. 
2015). Therefore, the disparity between dark gas-rich galax¬ 
ies (without clear stellar counterparts at present, e.g. 
Tollerud et al. 2015), and the paucity of gas detected 
near the stellar overdensities identified by K15 and B15 
(Westmeier et al. 2015), can be traced to the highy ec¬ 
centric/energetic trajetory of the Magellanic Clouds/Group 
which penetrates some of the densest regions of the Galactic 
hot halo. 

By accounting for the possible influence of the Mag¬ 
ellanic Stream, the primary uncertainty pertaining to our 
ram pressure model lies in the gas stripping rate. Nonethe¬ 
less, our second key prediction lies in observable signatures 
of this mechanism, in the form of compact HVCs lying in 
the wake of the Magellanic Stream and distinguishable by 
their low metallicity. Accordingly, the UFDs objects are ex¬ 
pected to be dominated by an old stellar population with 
a negligible specific star formation since commencement of 
their recent accretion between 1 and 5 Gyr ago. 
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